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ABSTRACT 


<  i 


,> Mathematical  models  used  in  design  of  the  XM204,*Howitzer 
(AD  and  ED  Prototypes)  are  described, in  this  report,-'  The  physical 
basis  for  the  mathematical  representation  is  presented  along  with 
the  derivation  of  the  required  equations.  While  these  models  have 
been  generalized  to  allow  their  use  in  other  weapon  design  situations, 
some  modification  will  be  necessary  to  include  features  not 
specifically  considered.  Systems  of  equations  which  will  provide 
for  the  definition  of  required  control  functions  as  well  as 
the  prediction  of  recoil  mechanism  functioning  and  weapon  motions 
are  summarized. 


FOREWORD 


This  project  was  authorized  under  AMCMS  Code  553G.12.42710.01. 
The  work  was  funded  under  DA  Project  No.  1-W-5-63608-D-376 . 

Suggestions  and  requests  by  the  members  of  the  XM204  design 
team  served  as  the  basis  for  the  mathematical  representations 
of  the  weapon  and  their  support  is  hereby  acknowledged. 
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INTRODUCTION 


The  recoil  mechanism  of  an  artillery  weapon  provides  a 
controlled  resistance  to  some  allowable  motion  of  the  recoiling 
parts.  Consequently,  the  major  portion  of  the  weapon  structure 
is  protected  from  the  total  breech  force  and  the  ground  reactions 
are  limited.  In  a  conventional  recoil  cycle,  motion  is  initiated 
by  application  of  the  breech  force.  The  resisting  force  may  be 
provided  by  deflection  of  a  spring,  by  fluid  flow  through 
an  orifice  or  by  both  of  these.  In  the  soft-recoil  cycle,  the 
recoiling  parts  are  first  accelerated  in  the  direction  of  projectile 
travel.  Application  of  the  breech  force  (i.e.,  by  ignition  of 
the  propellant  primer)  is  delayed  until  a  predetermined  velocity 
has  been  attained.  The  breech  force  then  reverses  the  motion  of 
the  recoiling  mass  and  the  recoiling  parts  are  returned  to  their 
initial  (battery  or  latch)  position.  A  comparison  of  conventional 
and  soft-recoil  cycles  is  shown  in  Figure  1. 
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FIGURE  1 


Comparison  of  Recoil  Cycles 


Development  of  a  reliable  soft-recoil  mechanism  for  modem 
artillery  began  in  1957  with  the  modification  of  an  M101  Howitzer, 
105om.  While  evaluation  tests  of  this  modified  weapon  demonstrated 
the  feasibility  and  value  of  a  soft-recoil  system,  mechanism 
reliability  was  unacceptable.  (Reference  1),  In  March  1964, 
design  and  fabrication  of  an  experimental  firing  fixture  having  an 
improved  recoil  mechanism  was  initiated.  Extensive  firing  tests  of 
this  fixture  were  conducted  to  confirm  feasibility  of  this  weapon 
concept,  to  determine  accuracy  and  durability  characteristics,  and 
to  identify  and  examine  functional  problems.  (References  2  &  3) , 

In  1968,  funding  was  provided  for  design  and  manufacture  of  a 
preprototype  howitzer  incorporating  a  soft-recoil  machanism. 

This  weapon  was  designated  as  the  Howitzer,  Light,  Towed:  105mm, 

Soft  Recoil,  XM204.  The  major  components  were  the  XM205  Cannon, 
the  XM46  Recoil  Mechanism,  and  the  XM44  Carriage.  Design  parameters 
for  the  first  prototype  were  based  on  firing  of  the  XM606  (28.5  pounds) 
projectile  with  the  M85  charge  (Zones  3  through  8)  and  the  standard 
Ml  (33  pound)  projectile  with  the  M67  charge  (Zones  1  through  7). 

Safety  Certification  Tests  were  conducted  at  Aberdeen  Proving 
Ground  in  early  1970.  Military  Potential  Tests,  performed  at 
Fort  Sill,  Oklahoma,  were  successfully  completed  in  December  1970. 

Both  standard  I05im  howitzers  (M101A1  and  M102)  were  used  for 
comparison  during  these  tests.  Stability,  accuracy,  and  human 
engineering  characteristics  of  the  XH204  were  favorably  commented 
upon  by  the  user.  In  all,  2,269  rounds  were  fired  from  this  MPT 
weapon  with  413  at  the  maximum  impulse  level.  (Zone  8). 

(References  <4,  5,  and  6). 

Mathematical  models  of  the  weapon  were  developed  to  aid  the  design 
engineers  in  establishing  and  evaluating  physical  configuration, 
structural  integrity  and  functional  controls.  (References  7,  8  &  9). 
These  models  were  used  to  verify  designs  and  to  systematically  study 
the  effects  of  parametric  variation  before  selection  of  specific  values 

In  the  early  part  of  1972,  the  design  of  the  XM205  Cannon  was 
changed  to  permit  test  firing  of  the  XM20G  supercharge.  This  longer 
and  heavier  tube  was  installed  on  the  original  prototype.  Concurrently 
a  new  velocity  sensor,  an  improved  suspension  system,  ring  spring 
assemblies  in  the  elevation  struts,  and  a  relocated  firing  base  were 
incorporated.  (Reference  10).  With  fabrication  and  installation 
of  a  new  cradle,  this  modified  weapon  became  ar  Advanced  Development 
Prototype. 
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Design  of  the  Engineering  Development  Prototypes  that 
will  be  used  for  DT  T370T  XE  has  been  initiated.  As  a  result 
of  the  changes  in  the  cannon  and  the  increase  in  ammunition  impulse 
levels,  new  recoil  and  counterrecoil  orifice  designs  were  required 
and  the  effect  of  carriage  flexibility  had  to  be  reassessed. 
Consequently,  the  original  mathematical  models  were  modified 
to  more  closely  represent  the  current  design  configuration.  At 
the  same  time,  an  attempt  was  made  to  generalize  the  form  so  that 
minor  alterations  would  allow  for  their  use  in  the  analysis  of 
future  soft-recoil  weapon  concepts. 


CONCLUSIONS  AND  RECOMMENDATIONS 


While  the  models  presented  in  this  report  were  being  developed, 
the  primary  considerations  were  choice  of  significant  motions  and 
physical  characteristics  of  the  weapon.  Specific  design  features 
and  functional  characteristics  peculiar  to  the  XM204  Howitzer 
have  been  included  to  adapt  the  models  to  the  requirements  of  the 
design  team.  These  models  ate  intended  to  provide  a  reasonable 
representation  of  normal  firing  cycles  based  on  firing  of  the  standard 
zoned  charges  and  of  abnormal  cycles  resulting  from  a  cook-off 
or  from  a  misfire. 

Test  plans  include  the  collection  of  data  from  the  firing  of 
the  Advanced  Development  Prototype  under  various  conditions. 

Comparison  of  test  data  and  predictions  obtained  from  the  mathematical 
models  will  provide  a  basis  for  model  validation.  The  models  should 
then  be  used  to  identify  critical  design  parameters,  to  define 
the  effect  of  parameter  variations,  and  to  establish  required  values 
for  specific  parameters  for  use  in  design  of  the  Engineering 
Development  Prototypes. 
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Figure  6 

Fluid  Flow  Diagram 


RIGID  BODY  MODEL  OF  A  SOFT  RECOIL  MECHANISM 


The  schematic  diagram  of  the  soft  recoil  mechanism  (Figure  2) 
and  the  free  body  diagrams  (Figures  2,  3,  and  4)  are  the  basis 
for  the  mathematical  model  of  recoil  motion  used  to  predict 
fluid  flow  and  its  control.  Conventionally,  the  breech  force  is 
assumed  to  cause  positive  acceleration.  Therefore,  the  displacements 
"x",  "y"  end  "x  -  y"  increase  in  magnitude  for  the  force  system  shown. 
With  the  recoiling  parts  in  their  initial  or  battery  position, 
the  variables  "x"  and  "y"  are  defined  as  zero.  (Note  that 
"x  y  0"  in  Figure  2  -  i.e.,  rearward  displacement).  Symbols 
used  in  ueriving  the  equations  used  to  predict  system  motion  and 
functioning  a defined  in  the  Symbol  Table  (Page  80), 

The  spear  buffer  shown  in  Figure  2  is  required  to  protect  the 
system  from  the  overload  which  would  be  caused  by  the  firing  of 
the  maximum  impulse  charge  before  imparting  a  forward  motion  to  the 
recoiling  parts  (as  in  the  case  of  a  cook-off).  This  spear  buffer 
restricts  fluid  flow  between  the  pressures  P  and  P^  during  part 
of  the  cycle.  As  illustrated  schematically  in  Figure  5,  this 
restriction  exists  while 


with  "xe"  being  the  value  of  "x"  for  which  a.^  becomes  equal  to  . 

For  any  "x",  the  flow  diagram  is  shown  in  Figure  6  with  flow 
from  P^  to  when  "x"  is  increasing  (defined  as  recoil).  The  direction 
of  flow  is  reversed  when  "x"  is  decreasing  (during  run-up  and  during 
counterrecoil).  Values  for  the  orifice  areas  (a/s)  will  vary  with 
direction  of  fluid  flow,  position  of  the  spear  buffer  relative  to 
the  piston,  and  functioning  of  the  velocity  sensor. 


The  pressure  drop  across  the  "ith"  orifice  is  defined  as 

o/M1 

(1)  g(v.)  -  —  I  —  I  sgn(Vi) 

Then, 


(2) 

P1  - 

P2  * 

8^) 

(3) 

F2  " 

P3  - 

g(v2) 

(4) 

P3  - 

P4  ■ 

gW3) 

7 


wile  re 


O’-  fluid  density 
3  —  acceleration  due  to  gravity 

*  discharge  coefficient  for  "ith"  orifice 


sgn(v^) 


\  ”  fluid  velocity  through  "ith"  orifice 

)  *  algebraic  sign  of 


With  fluid  flow  restricted  by  the  spear  buffer  (that  is. 
for  "x  £  x  ’ 


NA^x  *  Na2vi 


^1 

71  ax  x 


(6)  v 


Nalvl 

+  nA2*  * 

Na2v2 

NA  * 

+  na2x  « 

Na2v2 

A1  +  A2  . 

v2  " 

x 

Na2v2 

=  NA,x  + 

-a 

a3V3 

a3V3 

=  NAj  x 

NA^x  -  NA^x  =*  NAjjX 

(7)  v. 


~r -  x 


Nov 


a3v3  -  A^(x  -  y) 

NV  “  V  "  V 


(8) 


-  nar 


(9) 


•  • 

y 


'  nar 


(10) 


y  * 


A4-  NAR 


x 


These  relations  hold  as  long  as  the  system  remains  completely  filled 
with  fluid.  (Monitoring  of  computed  pressure  values  to  ensure  that 
they  remain  positive  will  demonstrate  the  existence  of  an  oil 
filled  system) . 


Rearranging  Equations  2,  3,  and  4 

(11)  P3  -  P4  +  g(v3) 

(12)  P2  =*  P3  +  g(v2)  =*  P4  +  g(v2)  +  g(v3) 

(13)  pi  “  p2  +  8(v1)  ■  P4  +  g(v^  +  g(v2)  +  g(v-j) 

From  the  free  body  diagrams  (Figures  2  and  3) ,  with  "x  »x,"  (spear 
buffer  restricting  fluid  flow) 

(14)  x  «  B(t)  +  WR  sin  ^  -  NFp  sgn(x) 

-  Ffp  sgn  (x  -  y)  -  Fg  sgn(x) 

~  Vs  +  A4P4  +  NA3P3  -  NA2?2  *  NA1P1 
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and 


(13)  M?  y  -  W  sin  J  +  +  Pf  sgn(i  -  y)  - 


Noting  that  sgn(x)  *  sgn(t  -  y)  and  substituting  for  y  (Equation  9), 
Equations  14  and  15  are  rewritten  as 


(14.1)  MR  x  -  S(t)  +  WR  sin  -  (NFp  +  Fg  +  Ffp)  sgn(x) 

'  A/n  +  V*  +  NA3P3  ”  NA2P2  "  NA1P1 


(15.1) 


A,  -  NA_*| 

Mp  - A^H  *  “  WP  Sln  J  +  Vn  +  Ffp  S*n(i)  ‘  AAP4 


By  adding  Equations  14,1  and  15.1,  one  obtains 


C16)  |“* +  |— Mpj  *  ■  “"=>  +  (»* +  y  sl»7 

-  (NFp  +  F  )  sgn(x) 

+  NA3P3  -  NA^  -  NA1P1 


By  substituting  Equations  11,  12  and  13,  Equation  16  may  be  rewritten  as 

["R +  (■**“)  M»]  v " B<t)  +  Wr  +  V  sinf 

-  (NF  +  F  )  sgn(x)  +  NA,P. 

p  g  3  4 

+  NA3  g(v3)  -  NA^P ^  -  NA2  g(v2) 

-  NA2  g(v3)  -  NA^  -  NA3  gCvj^) 

-  NA^  g(v2)  -  NA1  g(v3> 


-bi’lAsfe.-. 


Collecting  terms 


k +  (-4^j  \  *  ■  ,(t>  +  *  V  *ln  7 


-  (HFp  +  Fr)  agn(i) 


-  ^  -  A^)  P4 

-  g(vx)  -  g(v2) 

-  NCAj.  +  ^  -  A3)  g(v3) 


Solving  Equation  15.1  for  P^, 


(15-2>  p«  •  4;"p  •**?*  4:  %  a8n(*) 


+  - -  P 

ku  *N 


a4  -  nar| 


By  substitution,  (noting  that  Ap  »  “S.  +  A2  ~  A3^ 


k + 1~ a7~ ]  Mpj v  ■ B(t)  +  ("r  +  V  *1”  7 


NA_ 

-  (NT  +  F  )  sgn(i)  -  —■  W_ 


gCv^  -  +  Aj)  g(v2)  -  HAg  g(v3> 


or,  rearranging  terus 


(17) 


x  -  B (t)  + 


[- 


'h  -  NM 

A4 


NAs  .  naran 

-  (NF  +  F  +  -jS  F  )  sgn(x) - 

P  g  A,  £p  A, 


"pi  ■1"  J 


NAj^  g(v1)  -  N(A^  +  A2)  g(v2)  “  ^  g(v3) 


If  the  initial  value  if  V  ,  the  gas  volume  for  any  displacement  "x" 
is  written  as 


v»  ‘  vo  -  S  <*  -  y> '  vo  -  *N  r*  » 

4 

Then,  assuming  adiabatic  gas  laws,  the  gas  pressure  for  any  displacement 
"x"  is  determined  from 


P  V 
N  N 


(18) 


From  the  free,  body  diagram  (Figure  4),  the  force  cn  the  recoil  rod  is 
given  by 


(19)  R  -  NA.P.+  NA.P,  -  wA_F.  NF  sgn(x) 
IX  2  Z  3  3  p 


When  fluid  flow  is  not  restricted  by  the  spear  buffer,  (i.e. 
for  <  x  "),  the  preceding  analysis  a\ist  be  modified  by  noting 
the  force  on  the  end  of  the  spear  buffer  (Figure  2)  will  be  NA-P^ 
rather  than  NA2?„  and  that  flow  from  to  P^  will  be  defined  By 

N(A^  +  A^  x  *  Naivj 


(5.1)  vx 


Al  +  A2 

-  x 

al 


that 
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Therefore,  collecting  terns  and  noting  that  A^  •  A^  +  -  A^ 


[s-M 


B(t)  +  (WR  +  W )  sin J 


(NFp  +  F  )  sgn(i) 


-  NA^  -  N(A1  +  Aj)  g(vL) 


-  N(A1  +  A^  g(v2)  -  g(v3) 


Again,  substituting  for  and  rearranging. 


h*!^)  . |vM 


2 


(17.1) 


% 


-  (NFP  +  Fg  +  Af  V  8*“(X) 


nVn 


N(A3  +  A2)  [g(v1>  +  g(v2>] 


"  NAR  g(v3) 


By  defining 


Equations  S  and  5.1  may  be  written  as 


-  +  HAj)  g(Vl) 

-  N(A^  +  A^  g(v2> 


-  NAr  g(v3) 

Then,  by  using  Equations  5.2  and  18.2  along  with  Equation  20, 
a  single  system  of  equations  defines  the  entire  cycle  as  long  as 
the  mechanism  remains  filled  with  fluid. 

Negative  pressure  values  obtained  during  the  solution  of  the 
preceding  system  of  equations  will  denote  that  the  system  is  no  longer 
completely  filled  with  fluid.  If  this  occurs,  a  new  set  of 
equations  is  required.  This  condition  is  anticipated  since  fluid 
flow  through  the  velocity  sensor  (orifice  area  a  )  is  sharply 
restricted  after  Initiation  of  firing  to  preventJany  Increase  in  velocity 
during  an  ignition  delay  period.  To  approximate  system  motion  after 
a  negative  pressure  has  been  computed,  assume  that 


and  Equations  14  and  15  must  be  rewritten  as 


(21)  x  -  B(t)  +  WR  sin  -  (NFp  +  F  )  sin(i) 

-  Ffp  *»gn(x  -  y)  -  /^PN  -  A4  g(v3) 

(22)  Mp  y  *  Wp  sin  ^  +  Ffp  sgn(i  -  y) 

+  A4  g(v3) 


with 


a3v3  “  \  ~  y> 


or 


(23)  v. 


(x  -  y) 


and 


(24) 


P„  *  P 
N  o 


Vo  “  (X  “  y> 


These  equations  will  be  considered  as  definitions  of  system 
motion  as  long  as 

NA^  x  >  A4  (x  -  y) 

During  this  period 

(25)  R  -  NFp  sgn(x) 
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The  governing  differential  equations  of  motion  are  solved 
by  standard  numerical  methods,  with  the  system  controls  of  the 
actual  mechanism  being  simulated  by  logic  decisions  of  the  computer. 
These  logic  decisions  are  defined  in  the  following  manner. 


1.  \  cycle  is  initiated  by  the  setting  of 

appropriate  initial  conditions  (e.g.,  x  »  0 
and  x  *  0  at  t  ■  0). 


2.  Firing  is  initiated  on  the  basis  of  either 
velocity  (Is  "x”  i  some  specified  value?) 
or  displacement  (Is  "x"  some  specified 
value? ) 

3.  Ignition  delay  is  simulated  by  a  variation 

of  the  rime  lapse  between  initiation  of  firing 
and  a^v^-i-cation  of  the  breech  force. 

4.  The  orifice  area  a^  may  have  either  of  two 
values. 


a 


3 


a^bf  before  firing  if  "x  <  Q" 
<  a3af  after  firing  if  "x  <  0" 


3b  f 


whenever  "x 


0" 


5.  The  orifice  3rea  ^'aa  a  single  value  for  this 
model. 

6.  The  orifice  area  a^  varies  with  both  position  and 
direction  of  fluidiflow,  since  it  is  dependent 

on  restriction  of  fluid  flow  between  P  and  ?2  by 
the  spear  buffer  and  by  functioning  or  the  check 
valve  in  the  spear  buffer.  Letting 


*  flow  area  through  check  valve 

g 

leak  «  annular  orifice  area  resulting  from 
the  necessary  clearance  between  the 
piston  head  and  the  spear  buffer. 

a  *  orifice  area  which  is  dependent  on 

position  of  the  spear  buffer. 
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The  orifice  area  nay  be  defined  os 


“leak 

if 

"x  <  x  '* 
e 

a,  , 

+  a  if 

"x  >  X  >  X  " 

leak 

i 

V 

v  e 

and  "x  >  0" 

“leak 

if 

"x  >  xv”  and  "x  > 

0" 

“leak 

if 

"x  >  x  "  and  "x  < 
6 

0" 

”“x- 

Aj"  for  "x 

<  xe"  the  definition 

slaplifies  to 

“leak 

unless 

"x  >  x  >  x  "  and 
v  e 

o 

A 

“leak 

+  av  when 

"x  >  x  >  x  "  and  "i  >  0" 

v  e 

7.  Define 


"x  <  xe" 


"x  >  x  " 
—  e 


8.  Choose  governing  systea  of  equations  on  the  basis  of 


P3  >  0  or  P3  <  0 


The  rigid  body  aodel  of  the  soft  recoil  nechaniea  is  suseaarized 
in  Table  I. 


18 


TABLE  I 


RIGID  BODY  MODEL  OP  SOFT  RECOIL  MECHANISM 


NA 

-  (NF  +  Fb  +  ~  F.  )  »gn(x) 
P  8  \  tp 


-  N(A1  +  HA^  g^)  -  N(A^  +  g(v2> 
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TABLE  I  (coat'd) 

(6) 

Aj^  +  HA, 

• 

V,  - 

2  *2 

X 

(7) 

*1  +  A2 

• 

v,  - 

X 

(10) 

3  aj 

A,  -  NA* 
y  '  *» 

X 

(15.2) 

P4  ‘  \  Sl” 

7  *  Tk  Ffp  S8n<;> 

♦  — 

p  .  i  p*  -jsJ 

N  A4  I  A4 

A4 

(ID 

P3  ’  P4  +  «<V3> 

L  J 

(12) 

P2  -  P3  +  g(v2) 

(13) 

P1  -  p2  +  g(vx) 

(IS) 

R  -  NA1p1  +  NA2P2  -  NA3p3  +  NFp  sgn(i) 

X  >  A,  (x  -  y) 


B(t)  +  W  sin  /  -  (NF  +  F  )  sgn(x) 
K  r  P  8 

-  Ffp  agn(x  -  y)  -  A^  -  A4  g(v3> 

H„  •1“  7  *  Vs  +  rf ,  •«»<*  -  *> 


+  A,  g(v,) 


TABLE  I  (cont’d) 
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DESIGN  OF  RECOIL  AND  COUNTEREBCOIL  CONTROL  ORIFICE 
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with  "x  <  0"  (during  counterrecoil) ,  this  nay  be  written 

(after  assuming  c^  c2  **  C3  ”  48 


M  x 
eff 


% 

(NF  +•  F  +  -~ 

PS  A, 


-  F  ) 
4  fP 


NVn 


<t  iv 

23c2  [  ^ 


N(AX  +  A^  (NA^) 
2  +  .  2 


As  indicated  in  this  equation,  the  gas  pressure  "P  "  produces 
the  negative  acceleration  while  effective  friction  and  Hydraulic 
throttling  produce  the  positive  accelerations.  The  velocity  of 
the  recoiling  parts  will  be  constant  if  the  acceleration  is  equal 
to  zero.  That  is,  if  the  following  relation  holds 


N(A1  +  A^ 


(NA^T 

2 


nar  nVn 

+  CNF  +  F  +  F.  )  -  — 


To  ensure  control,  the  friction  forces  can  be  neglected  and,  by 
proper  sizing  of  the  orifice  areas  (a.),  any  one  orifice  can  be  made 
to  limit  the  counterrecoil  velocity.  A  specified  value  for  the 
minimum  counterrecoil  velocity,  (x  .  )  can  be  used  to  define  a  value 
for  f"  through  the  following  ?eSation. 


OV' 


nVn 


PN  ~ 


N(A1  +  KA 

2  "  ~ 


To  limit  the  terminal  counterrecoil  velocity  (x^) ,  the  value  of 
"a^"  can  be  defined  by  the  relation 


,2  K*»  A 


Actually,  the  minimum  value  for  a^  Is  determined  by  the 
necessary  clearance  between  the  spear  buffer  and  the  recoil  piston 
head.  Then,  at  latch  position, 

al  “  aleak. 

and,  by  neglecting  terms  that  include 


a2  811(1  a3af 


DEFINITION  OF  RECOIL  AND  COUNTERRECOIL  CONTROL  FUNCTIONS 


Linear  motion  (x)  of  the  recoiling  mass  (M^^)  is  defined  by  a 
differential  equation  having  the  fora: 


(40)  Meff  x  -  A(t)  -  D(t) 


where  (Figure  7) : 


A(t)  «  Summation  of  forces  causing  positive  acceleration 


D(t)  -  Summation  of  forces  causing  negative  acceleration 
and  both  may  be  considered  as  functions  of  time. 


FIGURE  7 

Free  Body  Diagram  of  Recoiling  Mass 


Since  the  firing  cycle  begins  with  the  recoiling  parts  at  rest,  the 
initial  conditions  are: 

x  *  0  and  x  *  0  when  t  *  0 

If  a  definition  of  A(t)  and  D(t)  is  assumed,  a  graphical  interpretation 
of  Equation  40  at  any  instant  (t  M  t.)  is  given  by  the  following 
equations.  (Reference  8,  Pages  6  ana  7). 


Moment  Area  Equation* 


(41) 


(42) 


t  t± 

Meff  *t  -  t,  "/  A<c>  dt  -J  D(t>  dt 
1  0  0 

\h  \  -  H  -[‘i  *(c)  dt 


ci 

-  [ti  -  D(t>] j  D(t) 


dt 


where 


(43) 


f: 


A(t)  dt  ■  Area  under  curve  of  A(t)  froa  t  •  0  to  t  ■  t 


rci 

I  D(t)  dt  -  Area  under  curve  of  D(t)  froa  t  •  0  to  t  •  t^ 


(44) 

J  D(t) 

o 

(45)  | 

>i  -  *(')] 

(46) 

[tt  -  D(t>] 

A(t)  dt  «  Moment  of  area  under  A(t)  around  t  -  t. 


(46)  -  D(t)]  J  D(t)  dt  -  Moment  of  area  under  D(t)  around  t  »  t^ 


If  A(t)  la  prescribed  forcing  function  (i.e. ,  a  breech  force 
plus  an  effective  weight  component) ,  D(t)  will  be  a  control  function  by 
which  some  specified  motion  of  the  mass  H  „  ia  produced.  Equations  41  and  42 
can  be  used  to  obtain  the  required  definition  of  the  control  function  C(t). 

For  a  conventional  recoil  cycle,  the  breech  force  [B(t)]  is  applied 
at  t  ■  0  and  recoil  displacement  must  be  limited.  Then,  let 


D  •  Value  of  D(t)  at  beginning  of  recoil  (t  ■  0) 
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D  •  Value  of  D(t)  at  end  of  recoil  (t  *  t  ) 
e  r 

1  ■  Allowable  recoil  displacement  (x  ■  l  at  t  ■  tf) 

I  -  Area  under  B(c)  (Impulse) 

a  -  Location  of  centroid  of  area  under  B(t) 

W  r,  •«  Effective  weight  component 
et  r 


A  minimum  peak  value  for  D(t)  can  be  obtained  if  the  function  is 
held  constant  over  the  complete  recoil  stroke.  However,  since  an 
instantaneous  increase  in  D(t)  is  impossible,  a  reasonable  rise  and 
fall  time  mus*.  be  allowed.  Let 


»  Specified  rise 
Dr  ■  Constant  level 


and  fall  time  for  D(t) 
of  decelerating  force 


The  assumed  force  system  is  illustrated  graphically  in  Figure  8. 


FIGURE  8 

Assumed  Force  System  for  a 
Conventional  Recoil  Cycle 


28 


Equation*  41  and  42  aay  be  aolved  simultaneously  for  value* 
of  the  unknown*  dr  and  t^. 


This  defines  the  control  function  used  to  design  the  spear 
buffer  which  protects  the  systea  from  daaage  If  a  atzlaua 
lapulse  charge  is  fired  from  latch  (i.e. ,  cook-off). 

For  an  ideal  soft-recoil  cycle,  the  assuaed  force  systea  Is  that 
shown  in  Figure  9. 


FIGURE  9 

Aasuaed  Force  Systea  for  an 
Ideal  Soft  Recoil  Cycle 
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D(c)  -  D 


(Aasuaed  constant) 


x  *  0  and  x  •  0  at  t  ■  0  (Battery) 

x  •  0  end  x  ■  0  at  t  ■  tR  (Return  to  battery) 

x  -  i  and  x  ■  x^  at  t  ■  t^  (Firing) 

the  moment  area  equations  (Equations  41  and  42)  can  be  solved 

D  -  Required  drive  force 

-  Cycle  tine 

xf  -  Firing  velocity 


(2)  returned  to  a  point  slightly  to  the  rear  of  latch 
position  by  application  of  the  breach  force;  and 


(3)  returned  to  latch  by  tha  ga a  force. 


During  a  maximum  overload  cycle  (cauaed  by  firing  from  the  latch 
position,  the  recoiling  parts  are 

(1)  driven  fearvard  by  application  of  the  breech 
force; 

(2)  brought  to  a  atop  by  hydraulic  throttling  of  oil 
between  the  apear  buffer  and  the  recoil  piston, 
and, 

(3)  returned  to  battery  (latch  poaitlon)  by  expansion 
of  the  comprensed  gas  in  the  recuperator. 

In  both  cases,  counterrecoil  control  must  be  provided  to  prevent  impact 
against  the  latch  mechanism. 

For  the  soft-recoil  mechanism  (shown  schematically  in  Figure  2) , 
the  minimum  counterrecoil  velocity  la  limited  by  restriction  of  the  fluid 
flow  from  the  recuperator  through-the  orifice  area  a^.  (After  firing, 
a^  *  a3af  48  defined  by  Equation  37  when  "x  <  0.")  To  protect 
the  latch  mechanism,  final  counterrecoil  control  is  produced  by  use 
of  the  spear  buffer  to  restrict  the  oil  flow  between  P2  and  P^.  Use 
of  the  spear  buffer  for  both  recoil  and  counterrecoil  control  1b 
made  possible  by  the  check  valve  which  is  effective  while  "x^  >  x  >  0." 
During  recoil,  and  the  valve  is  opened  to  allow  flmv  through 

the  valve  orifice  (ay)  as  well  as  around  the  spear  buffer  (ax  +  • 

During  counterrecoil,  with  P2  >  and  the  valve  closed,  flow 
between  P2  and  is  restricted  to  the  annular  orifice  around  the 
apear  buffer  (ax  +  aieak^  *  ^See  Pa*e 

For  an  acceptable  counterrecoil  control  function,  Che  counterrecoil 
velocity  muat  be  changed  from  x^n  at  x  »  xv  to  x^  at  x  ■  Xj  by  throttling 
of  fluid  flow  between  P2  and  P^.  As  a  result  of  tha  restricted 
flow  between  the  recuperator  and  recoil  cylinders  through  the  orifice  a^  ^ 
(thereby  limiting  the  magnitude  of  x^a) ,  the  accelerating  force  is 
approximately  zero  when  control  of  fluid  flow  through  a^  is  initiated. 
Since  the  terminal  velocity  is  to  be  held  constant  from  x  •  x^  to 
x  ■  0,  the  accelerating  force  will  again  to  zero  when  x  ■  Xj. 
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By  substituting  Equations  2,  3,  and  4  in  Equation  17,  w*  obtain  the 
following  equation  defining  counterrecoil  notion  at  0*  Q.2. 

[i.e.,  B(t)  ■  0,  x  0  and  sin  ^  *  0] 


(52) 


2 

)  ">] 


(N?p  +  Pg 


fp 


) 


NAj (P^  -  P2) 


-  N(A1  +  Aj) (?2  -  P3>  -  NAs(P3  -  P4) 

with  the  significant  pressure  drop  across  the  orifice  a^. 


p  =»  p,  fts  p  a.  p„ 

rN  4  3  2 


and  the  above  equation  becomes 

<53>  js  +  +  ^ 

"Vs 

-“iP  PN  -  »VP1  ’  V 


V 


If  P^  ■  0,  the  accelerating  force  is  given  by 

NAR 

D  -  (N?  +  F  ~  F,  ) 

P  8  “P 

*  PN  +  NVn 

nae 

D  -  (HF  +  F  +  ?,) 

P  8  A4  fp 


Vh 

+  N(Ai-^)  pn 
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With  counterrecoil  control  occurring  near  the  latch  position  where 
PN  **  p  ,  a  Uniting  design  value  for  "D"  is  given  by 


Based  on  the  preceding  analysis,  the  counterrecoil  control  function 
must  satisfy  the  following  constraints 


D(t)  "  0  when  x  •  *  and  x  ■  x  . 

v  min 

D(t)  5  D  for  xy  <  x  <  0 
D(t)  •  0  when  x  ■  x^  and  x  »  x- 


The  shape  of  the  control  function  may  be  chosen  in  some  essentially 
arbitrary  fashion  so  long  as  adequate  flow  control  can  be  maintained. 

One  such  choice  is  shown  (Figure  10). 


FIGURE  10 

Counterrecoil  Control  Function 


T 


Then,  from  "t  ■  0  to  t  •  y  ",  notion  is  defined  by 


2!) 

Meff  *  "  T 


CR 


where 


M 


eff 


nar 


Vn 


Dci  *  %  +  V  V  +  N(A1  “ -if  >  -o 


NOTE:  t '  is  assumed  zero  when  counterrecoil  control  begins. 


Then, 


2D 


Meff  X 


CR 


+  K, 


D„_  3 

CR  t 


Meff  X 


+  K  t  +  K2 


Since 


x  ■  x  and  x  ■  x  .  at  t  ■  0 
v  nun 


The  constants  ox  integration  are 


K,  “  M  x 
1  eff  min 


K,  ■  M  - ,  x 
2  eff  v 
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Therefore,  at  t  -  j  ,  x  «  X  and  x  -  X 


M  ;  °CR  T2 

X  *  “"t~~  “T  *  M  *  X  . 

r*  i  **  efr  mxn 


dcrt 


M*ff  X  “  '3T  'T  +  Meff  ^  I  +  Meff  \ 


D„  T  xT 

CR  .  Xmla  r 

f  -  4.  x 


While  the  function  D(t)  is  decreasing,  the  equation  of  motion  may  be 
writ  ten  as  7 


Meff  ‘x*  "  °CR  "  -T~  c 


x  *  X  and  x  *  X  at  t  =»  0 


NOTE;  Here,  "t"  has  been  assumed  zero  -?hen  the  control  function  begins 
to  decrease.  6 


Meff  *  “  °CR  C 


2D  2 

-a  il  +  r 

T  2  *1 


f_2  D  3 

H  .  x  “  D__  f  k  t  +  v 

eft  CR  2  T  3  K1  C  +  K2 


and  the  constants  of  integration  are 


K1  '  “aft  * 


*2  '  M.«  * 


As  originally  specifiad,  when  "D"  returns  to  zero 


The 


x  “  x^.  and  x 


n,  at  t  -  j  , 


M  i  » D  i  -  za  r 

eff  *T  UCR  2  4 

T 


+  M 


T 

CR 


eff  4M 


'eff 


+  M  ..  x 
eff  min 


Heff  “  °2  +  Meff  Xmin 


(35)  Vn  *  *T  " 


D  T 
CR 


2M 


aff 


Meff  24 


D  _3 

CR  T_ 

3T  8 


D  T 

+  y.  +  M  -  x 

4  2  eff  min 


D  T 
CR 

24 


+  M  x,  —  +  H  -  x 
eft  min  2  etf  v 


Meff  *T  “  DCRT  ^8  “  24  +  8  +  24] 


+M.-X,  T+M-.x 
eff  min  eff  v 
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N>|H 


Equations  55  and  56  may  be  solved  simultaneously  for  values  of 
"T"  and  "x^"  when  values  for  "DCR,  Me£f,  x^,  and  x^"  are  specified. 

By  substituting  Equation  55  in  Equation  56 , 


<S7)  °C»  2 

—  T  '  K.ff  *T  1  '  M.ff  <*v  -  *T>  '  0 


Equation  57  can  then  be  solved  for  "T"  and  "x  .  "  determined  by 
solving  Equation  55.  10111 
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CARRIAGE  RESPONSE  TO  FIRING  LOADS 

In  this  model,  three  degrees  of  freedom  are  considered  as  shown 
in  the  schematic  representation  (Figure  11).  The  mass  centers  shown 
are  defined  as: 

Mjj  -  Mass  of  nonelevating  portion  of  the  weapon 

Mg  -  Mass  of  elevating  (but  nonrecoiling)  parts 

Ma  »  Mass  of  recoiling  parts 

The  three  degrees  of  freedom  may  be  defined  in  terms  of  the  three 
time-dependent  variables: 

-  Defines  the  position  of  MA  with  respect  to  the  elevation 
trunnion 

9  -  Defines  relative  rotation  between  Mg  and  Mt,  around  the 

elevation  trunnion 

?  -  Defines  rotation  of  M^  around  an  axis  of  rotation  fixed 

in  the  ground 

On  the  assumption  that  all  motions  being  considered  are  in  the 
plane  of  elevation,  only  planar  coordinates  will  be  required.  Therefore, 
dimensions  are  defined  in  the  positive  sense  according  to  the  following 
coordinate  systems  as  shown  in  Figure  11. 

O-BC  :  A  coordinate  system  fixed  in  the  ground  with  its  origin 
at  the  point  of  rotation  for  the  mass 

0-YZ  :  A  coordinate  system  fixed  in  the  mass  MD  (and  rotating  with 

it)  having  its  origin  at  the  point  of  rotation  for  this  suss. 

T-15  :  A  coordinate  system  rotating  with  the  mass  and  having  its 

origin  at  the  elevation  trunnion.  At  the  time,  t  •  0,  this 
system  has  been  rotated  through  the  elevation  angle  (Y  ) 
from  its  original  position  parallel  to  the  0-YZ  coordinate 
system. 

From  Figure  11  ,  the  following  relationships  may  be  derived: 


ya 

-  yt  + 

cos  (y  + 

9)  -  ;A  sin  (Y  + 

9) 

ZA 

-  zT  + 

n  ^  sin  (Y  + 

9)  +  ?A  cos  (Y  + 

9) 
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Ba  -  Yt  cos  *  -  Zt  sin  $  +  nA  cos  (y  +  0  +  4)  -  CA  «in  <Y  +6+4) 

CA  *  Yt  sin  4  +  ZT  cos  4  +  nA  sin  (y  +  6  +  4)  +  CA  cos  (y  +  6  +  4) 

Similarity: 

YB  ■  Yt  +  ng  cos  (y  +  6)  -  Cg  3in  (y  +  6) 

Zg  -  Zj.  +  Hg  sin  (y  +  6)  +  Cg  cos  (y  +6) 

Bg  «  Y^  cos  4  “  Zj  sin  4  +  Hg  cos  (y  +  0  +  4)  —  Cg  sin  (y  +  0  +  4) 

Cg  »  Yr  sin  4  +  Zj.  cos  4  +  n0  sin  (Y  +  0  +  4)  +  Cg  cos  (Y  +  0  +  4) 

And, 

Bg  -  Yg  cos  4  -  Zp  sin  4 

Cg  -  Yg  sin  4  +  Zg  cos  4 

Noting  that  nA,  S  and  4  are  tine  dependant  variables,  differentiation 
with  respect  to  tine  gives: 

*  • 

BA  -  -  4  [Yp  sin  4  +  Zg  cos  4] 

-  (0  +  4)  (nA  sin  (y  +  0  +  4)  +  CA  cos  (y  +  0  +  :)] 

+  nA  cos  (y  +  0  +  4) 

•  • 

C/l  “4  lYg  cos  4  -  Zg  sin  4] 

+  (9  +  4)  [nA  cos  (y  +  9  +  4)  -  CA  sin  (y  +  9  +  4)] 


+ 

*nA 

sin  (y  +  9  +  4) 

• 

% 

m  — 

* 

sin 

4  +  Zj.  cos  4  ] 

- 

(0  ■ 

+  4)  [ng  sin  (Y  +  9  +  4) 

+  Cg  cos  (y 

+  0  +  0 )  ] 

• 

-  4  I^r 

cos 

4 

-  Zj,  sin  4  ] 

+  (0 

•f 

4) 

[  HgCOS  (Y  +  9  +  4)  - 

Cg  sin  (y 

+  6  +  4)] 

% 

• 

• 

Iyd 

sin 

4 

+  Z  D  coe  4j 

• 

<b 

• 

-  * 

Iyd 

cos 

4 

-  Z  D  sin  4] 
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Differentiating  again  with  respect  to  cine  gives: 


-  4>  [Yt  sin  4>  +  cos  4>] 

-  (0  +  $)  [  sin  (Y  +  S  +  +  ?A  cos  (Y  +  9  +  $)] 


+  "A  cos  (Y  +  8  +  4>) 

-  (9  +  $)  [2nA  sin  (Y  +  9  +  $)] 

9  2 

-  $>  [yt  cos  ♦  -  sin  4>] 

-  (9  +  40^[nA  cos  (y  +  9  +  $)-  CA  sin  (y  +  0  +  $) ] 

-  V  [Yt  cos  *  -  Zj  sin  $] 

+  (6  +  V)  [  nA  cos  +  6  +  $)  -  CA  sin  (y  +  9  +  $)] 

+  **A  sin  (y  +  9  +  <p) 

+  (0  +  40  [2nA  cos  (y  +  e  +  $)] 

.2 

-  i  [Yt  sin  <(>  +  Zj  cos  $] 

-  (9  +  {nA  sin  (Y  +  6  +  $)+  cos  (y  +  9  +  40] 

*  -  4  [Yt  sin  0  +  Z^.  cos  4] 

-  (6  +  4)  [nB  sin  (y  +  9  +  4)  +  Cg  cos  (Y  +  9  +  4)] 

•  2 

-  $  [Yt  cos  t  -  Zj  sin  4] 

-  (0  +  4)^  [0g  cos  (Y  +  3  +  ♦)  -  Cj  sin  (Y  +9  ->•  4)] 

-  4  [Yt  cos  t  +  Zj  sin  4] 

+  (0  +  4)  [nB  cos  (Y  +  0  +  ♦)  -  Cjj  sin  (Y  +  8  +  4)J 

*  2 

-  4  IYX  sin  4  +  Zj  cos  4] 

-  (0  +  i)2  [  0B  sin  (Y  +  8  +  4)  +  t_  cos  (Y  +  9  +  4)  ] 


cos  $] 


BD  * 


CD  * 


-  «  [Yd  sin  *  +  ZD 

.2  , 

-  *  [YD  cos  ^  -  ZD  sin  $] 

V  [YD  cos  $  -  Zjj  sin  <J>  ] 

*2 

-  £  IYd  sin  <t>  -f-  ZD  cos 


By  defining 


T  »  Kinetic  energy  of  the  system 
V  ■  Potential  energy  of  the  system 

F .  ,F.  ■  Generalized  forces  (torques)  causing  rotations  $  and  d 

9  □ 

the  Lagrange  equations  become: 

For  pitch  notion  (defined  by  $  ) : 


d_  u 

dt 


+  XL  .  F 


For  relative  rotation  (defined  by  8) : 


fir],  il  +  iV  . 

[39J  39  36 


If  XAl  Ig  and  In  are  mass  moments  of  inertia  of  each  mass  around 
its  mass  center,  the  Kinetic  energy  may  be  defined  in  terms  of  velocity 
components  (referenced  to  the  fixed  coordinate  system  0-BC)  and  rotational 
velocities  as: 

1  *2  .2, 

T  -  J  MA  (Ba  +  CA  ] 

1  .2  .2, 

+  2  MB  [*B  +  CB  1 

1  .2  .2, 

+  2  1»d  +  cd  ^ 

+  7  fXA  +  XB^  t ®  +  *1 


I  XD 
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After  terms  are  expanded  and  collected,  the  following  expression 
for  kinetic  energy  is  obtained: 

(60)  T  “  “7  LYt2  +  ZT2  +  nA2  +  ?A2  +  2<nAYx  +  r^Zy)  cos  (y  +  6) 

+  2<nAZx  -  caYT>  (y  +  6)] 

•  2  2  2  .2 

+  0  [nA  +  ;A  1  +  nA 

+  2*6  [  n  2  +  c  2  +  <n  Y_  +  C  2_)  cos  (Y  +  6) 

A  A  A  i  A  1 

+  <nAzT  -  cayt>  8ln  +  0)] 

-  2:  e  n.  -  2n  *  [t  -  Y  sin  (y  +  0)  +  Z  cos  (y  +  9)) 

A  A  A  A  i  1 

+  [Yt2  +  ZT2  +  n32  +  CB2  +  2  (nBYT  +  CBZT)  cos  (y  +  6) 


+  2  (n^  -  t3YT)  sin  <Y  +  9)] 


.  2  2  2 

+  9  [n_  +  C,  ] 


+  2*0  [nB2  +  CB2  +  (nBYT  +  CgZj)  cos  (y  +  0) 


+  (tib2j  -  CbYt)  sin  (y  +  0)] 


+  -r  ♦  iv +  zd  1  + 


lA  +  S  ,i 


.  2  h 

I*  +  9]  +  -s*  * 


Differentiating  Equation  60  with  respect  to  $  : 


3T 

w 


Ma[|  [Vt2  +  ZT2  +  nA2  +  5A2  +  2(nAYT  +  C^)  cos  (y  +  0) 

+  2(nAZT  -  CaYt)  sin  (y  +  6)] 

+  ®  [nA2  +  ;A 2  +  (nAYT  +  CaV  C0S  (Y  +  9) 

+  (nAZT  -  CaYt)  sin  (y  +  9)] 


-  n.  [;  -  Y  sin  (Y  +  e)  +  ZT  cos  (y  +  9)] 

A  A  A 


] 


hJ;  (Yt2  +  ZT2  4  nB2  »  tB2  4  2(nBYT  4-  co.  (Y  4  0) 

+  2(n  fa  -  CbYt)  sin  (Y  +  0)] 


0  [n2  +  CB  +  <nBYT  +  CBZT)  cos  (y  +  9) 


+  <Vt  “  W  sin  (Y  4 


4  Mj,  J[Yd2  4  Zd21  4  IIA  4  IB1(»  4  S]  4  h  i 


8)] 
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And  then  with  respect  to  time 


<61>  d£  [if]'  V  ["*  l\2  +  ZT2+  "a2  +  'a2  +  2<Vt  +  W  "•  ty  *  *’ 

+  2  (n^T  ~  ?aXt)  sin  + 

+  Mg  [Yt2  +  Zt2  +  aB2  +  ;B2  +  2(nBYT  +  cos  ^  +  °) 

+  2(nBZT  -  CByt)  sin  <Y  +  6>] 

Md  [yd2  +  ZD2]  +  [IA  +  iB  +  Id] J 
+  V  [ma  [nx2  +  Sx2  +  (nAYT  +  SAZT)  cos  (Y  +  9) 

+  (nAZT  "  5xV  8in  (Y  +  9)3 

+  Mg  [nB2  +  ?B2  +  (nBYT  +  ?BZT)  co8  (Y  +  0) 

+  (?bZt  -  5BYT)  -  Sin  <Y  +  9) ] 

+  t1* +  V  ] 

-n  Ma  [  CA  -  yt  sin  (  r  +  9)  +  ZT  cos  (Y  +  9)  ] 

-  2  *^MAf(nAYT+  8itl  (Y  +  6>  ‘  <nAZT~  SAV  co8  (Y  +  6)3 

+  Mb  [nBYT  +  c3ZT)  sin  (Y  +  9)  -  (nBZT  -  CgY.p  cos  (Y  +  9>] 

+  2®  MA  tnA  +  YT  C0S  ^  +  0)  +  ZT  Sin  (Y  +  9^3 

+  2  Ma  [nA  +  Yt  cos  (y  +  9)  +  ZT  sin  (Y  +  9)] 

-  92jtiA  [nAYT  +  ^ZJ  sin  (Y  +  9)  -  <nAzT  “  ^Y^  co8  <Y  +  6>  ] 

s 

+  Kb  [nBYT  +  sin  (Y  +  9)  -  (nBzx  -  5^  cos  (y  +  9)] 


Similarily,  differentiating  Equation  60  first  with  reapect  to  6 


■H*[i 


1  %2+  O 


+  *  (nA2  +  CA2  +  (Vt  +  fAZT)  CO“  <Y  +  e>  +  ^AZT  “  «ayt)  ^  +  <>)  1 

-  <A] 

+  ^  l^B2  *  S2' 

+  ♦  [nB2  +  CB2  +  (nBYT  +  ?BZT)  cos  <Y  +  e>  +  (nBZT  -  C^)  ain  <Y  +  0)] 


+  UA  +  ib3  [♦  +  e] 


And  then  with  respect  to  time, 


d  f  if!  eaf 

dT  [tfj"  *[“. 


A  ‘"A2  +  'a2  +  <Yt  +  W  +  »> 


+  (nAZT  -  CaYt)  ain  <Y  +  6)] 


+  **8  [nB2  +  S2  +  <Vt  +  W  C0#  <Y  +  9) 


+  (nBZT  -  tBYT)  .In  Or  +  8)1 


A  [IA  A  IB] 


+  V[ha  I8A2  +  'a2)  +  HE  l  V*  <b!  +  tl*  *  I,1]  +  \  1  -  W 

-  J»f«A  [(  "A«T  1  IAZT)  sin  (-r  +  8)  -  (nAz,r  -  4a(t)  cos  Or  +  e)j 

+  *B  C(n  BYT  +  CBZT)  8in  <Y  +  0)  ”  (  nBZT~  SV  008  <r  +  6)]] 

+  enA  [2Ma  nA]  +  nAJ  «A  {2nA  +  Yt  cob  (y  +  0)  +  z ain  (y  +  0)j 
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Next,  differentiating  Equation  60  with  respect  to  ♦ 


<63>  $  "  0 
and  with  respect  to  6 

(64)  — |I—  -  -  5»e^MA  I<nAYT  +  sin  (Y  +  8)  -  (n^  -  cos  (y  +  8)] 

+  Mb  [(nBYT  +  sin  (Y  +  8)  -  (nfiZT  -  C^)  cos  (y  +  8)] 

+  enA  [0] 

+  n,  J  M.  [Y_  cos  (Y  +  8)  +  Z_  sin  <y  +  8)] 

A  A  T  T 

-  ^2^Ma  [(nAYT  +  W  sin  (Y  +  8)  -  (nAZT  -  coe  (y  +  6)J 

« 

+  Mg  tCnBYT  +  CBZT)  sin  (Y  +  »)  -  (l^  -  C^j)  coe  (y  +  8)] 
The  potential  energy  in  the  syatea  aay  be  defined  as  the  aua  of 


Vy  -  Potential  energy  due  to  weight 

VG  ■  Potential  energy  stored  in  the  effective  spring  at  the  front  support  point. 


V 


EQ 


Potential  energy  stored  in  the  equilibratovrs 


Potential  energy  stored  in  the  elevation  struts 


Thac  is. 


+  V„  +  voo 
G  EQ 


Since  the  change  in  potential  energy  with  respect  to  each  of  the 
generalized  coordinates  is  the  important  factor  in  the  Lagrange  equation, 
any  reference  position  nay  be  used.  This  analysis  is  based  on  the 
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following  definition  for  the  reference  position: 


<j>  •  0 

0=0 


The  potential  energy  due  to  weight  is  given  by: 


Vw  “  CA  WA  +  CB  WB  +  CD  WD 


or 


(66)  V  »  [Y  sin  4>  +  Z_  cos  $  +  1,  sin  (Y  +  +♦)+?,  cos  (Y  +  9  +  $)]  if 

W  T  T  A  A  A 

+  [Y  sin  $  +  Z  cos  $  +  n  sin  (Y  +  9  +  <j>)  +C  cos  (y  +  0  +  <J>)  ]  W 
XT  d  fi  d 

+  [Yd  sin  <f>  +  ZQ  cos  $]  WD 


Then, 

(67) 


(68) 


differentiating  Equation  66  with  respect  to  $  and  with  respect  to  9, 

3V 

=  [Y_  cos  $  -  Z_  sin  <t>  +  n  cos  (Y  +  6  +  $)  -C  sin  (y  +  9  t  $)]  W 
cy  T  T  A  A  A 

+  [Y_  cos  $  -  Z,_  sin  $  +  n_  cos  '  9  +  c  '<  -  c  sin  (y  +  9  +  ♦)]  W 

T  i  3  3  o 

+  [Yd  COS  $  -  Z3  sin  ?>j  VID 

«V 

=  [h  cos  (y  +  9  +  4i)  -  ;  sin  (Y  +  9  )) 

+  [nB  cos  (Y  +  9  +  $)  -  cB  sin  (y  +  9  +  $)]  Wg 
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Lotting 


9tatic 


static 


3a  wa  +  bb  wb  +  bd  wd 


static 


The  potential  energy  can  be  written  as 


(69)  Vc  ■  i  KgL2  <♦  -  *  statlc)2 


Then,  differentiating  with  respect  to  41  and  2  gives 


V-  <♦  -  Static' 


The  equilibrators  and  elevating  struts  will  act  as  co-linear 
spring  elements  whose  potential  energies  change  with  their  change 
in  length.  To  establish  this  change  in  length,  the  following 
coordinates  may  be  defined: 


Y1  "  YT  +  nEL  C0S  <Y  +  e>  *  f'EL  sin  +  e) 


i 


Z  -  ZT  +  sin  (y  +  6)  +  cos  (y  -r  ?> 


Y  »  Y 

2  EL 


Z„  *  Z^r 
2  EL 


Then,  the  length  of  the  elevating  struts  is  given  by: 


-  V2  +  <Z1  -  Z2>2  1  1/2 


Substitution,  expansion,  and  collection  of  teras  allows  the  writing  of: 


LEL  t  Q  -  i  (Y  -  Y_  )  2  +  (Z_  -  Z„  )  2  +■  n  2+  2 

y  +  0  It  el  t  el  el  el 


*  2i(vt  -  V  ”il  +  <zi  -  ZEL>  W  +  "> 


2i«I  -  V  'jl-  <zr  '  ZEL>  "el1  *ia  ,r  +  s> 


By  defining  the  constants: 


(72)  ELI  -  <Yt  -  Yel)2  *  <ZT  -  ZEL)2  +  nEL2  +  ?EL2 

(73)  EL2  -  l(YT  -  YE1)  nEL  +  (ZT  -  ZEL)  ?EL] 

(74)  EL?  -  [<Yt  -  YgL)  -  <ZT  -  ZR)  nRi; 


The  length  may  be  written  as 


(75)  LEL  ^  +  „  -  [ELI  +  2EL2  cos  (y  +  0)  -  2EL3  ain  (Y  +  0)] 


1/2 
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A  spring  equilibrator  (compression  type)  in  shown 
schematically  in  Figure  13. 


FIGURE  13 

Spring  Equilibrator  -  Compression  Type 


Let 


PEQ  ■  Preload  in  equilibrator  spring  (i.e.  -  load 


■  Spring  rate  of  equilibrator  spring 


Then,  the  total  force  for  two  parallel  equilibrators 
is  given  by 


(76)  FEQy  +  g  -  2  [PEQ  -  K  (LEL,  +  e-  LELq)] 


With  the  reference  position  defined  by  9  *  0,  the  potential 

energy  stored  in  the  equilibrator  is  given  by 


V  -f 

EQ  J 


y  +  6 


feqy  +  6  d  (LEL  Y  +  e  ) 


-  j  2l'pEQ  -  keq  CLELy  +  6  -  LELj )  ]  d  (LEL  y  +  0) 


r  y  +  0 

2 J  [PEQ  +  Kjjq  LEI^  -  Kgq  (LEL  Y  +  a)]d  (LEL  Y  +  e  ) 

J  LEL.. 


LEL  .  .LEL 

f  Y  +  0  1  Y  +  0 

2  [(PEQ  +  Kgp  LEl^)  (LEL  Y  +6  )  -  KgQ  - *3 - J 


2  [(PEQ  +  Kgq  LEL^j )  (LEL 


-  LEL  ) 


EQ  2  2i 

2  (LELy  +  g  ~  LEL  Y  )J 


V 

8  EQ 

36 


-  2 [(PEQ  +  Kgq  LELg)  - 39 


Y  +  6 


°L£L  gl 

'  LELy  +  9  39  J 
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(77) 


V 

3  EQ 

38 


2  [PEQ 


-  LELq)] 


K__  (LEL 
EQ  Y 


3LEL 


yJJ 


36 


and,  since  V 


is  indpendent  of  the  variable  $ 


(78)  *  0 


The  elevation  struts  contain  preloaded  ring  spring  assemblies 
as  shown  in  Figure  14. 


FIGURE  14 

Schematic  Diagram  of  Elevation  Strut 
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Treating  an  increase  In  length  aa  a  p< 
spring  asaeably  and  defining  a  tension  lost 
one  to  depict  graphically  the  force  in  one 


Load  in  One  Strut 


H 

2 


P  r 

3 

2 


P 

3 


-  k(a  -  r)  +  P 

s 


PEL 

2 


Since  the  load  in  the  elevation  strut  must  be  defined  in  a  piecewise 
fashion,  the  potential  energy  will  also  be  defined  in  this  Banner. 

With  the  equilibrators  approximately  balancing  the  tipping  parts  when 
6  «  0,  the  potential  energy  in  the  elevation  struts  (both)  aay  be 

defined  as  twice  the  area  between  the  load  deflection  curve  and  the 
A  axle  from  A  »  0  to  A  The  curve  chosen  will  be  dependent  on 
whether  the  spring  assembly  is  being  compressed  or  releasing  its  energy. 

The  potential  energy  (V^)  nan  be  expressed  as  a  function  of  the 
strut  deflection  A  in  the  form 

C79)  VEL  "  14  “  r]2  K  +  2[*  '  r]  Ps  +  H 


where  the  constants  r,  K,  Ps  and  H  are  defined  in  a  piecewise  fashion. 
(See  Figure  15). 

Since  A  is  independent  of  <t> 


(80) 


-  0 


However , 


and  since 


,  3A  .  3A 
K  36  +  2ps  36 

l*  ~  +  Ps]f 


A  -  LELy  0  -  LELy 


The  equation 

3vel 

<81>  — 


[LEL 


Y  +  6 


lely 


r] 


+ 


3l:sl 


30 


+  0 


is  obtained 
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Recalling  that  (Equation  65) 


v-v+v+v  +  v 

W  G  EQ  EL 


The  following  relations  may  be  written 


34)  94*  9  $  3  $ 


3\  +  iVC  .  WEQ  „  WEL 
30  36  36  36 


Then,  making  the  appropriate  substitutions  in  the  above  equations,  using 
Equations  67,  70,  78  and  80 

(82)  —  ■  [Yt  cos  $  -  ZT  sin  *  +  nA  C03  (y  +  0  +  $)  -  ;A  sin  (y  +  6  +  <p>]  »A 

+  [Yt  cos  *  -  ZT  sin  *  +  nB  cos  (y  +  6  +  <fr)  -  ^  sin  (y  +  6  +4,)]  Wfi 

+  [ID  “S  *  -  ZD  Sl”  *  1  »„  +  KC  L2<*  -  »,t.tlc> 

Using  Equations  68,  71,  77,  and  81 


“  [nA  cos  (y  +  6  +  b)  -  CA  sin  (y  +  6  +  $)] 


+  [n*  cos  (y  +  6  i-  ♦)  -  U  sin  (y  +  6  +  b)  ]  wp 


♦  2  iPEQ-  Krq  (LELy+  9  -  LEL^)] 


+  2  [K(l,ELy  +  0  -  LELy  -  r)  +  Pg] 


which  may  be  written 
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By  the  resolution  of  the  breech  force  [B(t)]  into  a  force  along 
the  n  -  axis  end  a  couple  around  the  point  T  and  by  including  the 
dissipation  functions  representing  daaping  in  the  ground  spring 
and  in  the  elevation  trunnion,  the  generalized  forces  say  be 
taken  as  the  torques  around  the  appropriate  pivot  points  (0  and  T). 
Then, 

F  -  B(t)  r  -  Y  B(t)  sin  (y  +  e)  +  Z  B(t)  cos  (y  +  6)  -  C  9 
9  1  T  T  9 

F*  “  B(t)  ci  -  ce  e 

or 


(84)  F  -  Ic,  “  T  sin  (y  +  0)  +  Z  cos  (y  +  0)]  B(t)  -  c  $ 

9  1  T  r  f 

(85)  Fe  -  B(c>  -  e0  5 

Then,  since  the  Lagrange  Equation  for  the  pitch  motion  (defined  by  9)  is 


and  for  the  relative  rotation  of  masses  and  (defined  by  #)  is 


proper  substitution  enables  the  following  equations  to  be  written 
Equation  89  (See  Page  64  ) 

Equation  90  (See  Page  64  ) 
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Qj.  ■  ♦  [Yt  cos  (Y  +  9)  4-  Zt  sin  (Y  +  9)  + 

+9  CnA] 

-  2$9  CA  +29  nA  +  2  nA  % 

•  2  *  2 

-  8  4  -  $  [4  -  Y  sin  (Y  +  9)  +  Z  cos  (Y  +  9)] 

A  A  A  T  T 
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Mow , 


from  Figure  12(a)  (page  40) 

(86)  2Fn  -  Ma  \  -  R(t)  -  B(t)  +  +  f2  +  IQ1 

-  3F  -  W4  sin  (y  +  0  +  #) 

A 

(87)  ZTr  -  Ma  Q  -  S.  +  S,  -  W.  con  (y  +  0  +  $) 

^  A  r,  1  Z  A 

<88>  E!Ws  Center  '  E*  (V  +  ♦>  ‘  S2  ®  '  S>  '  S1  ® 

+  B(t)  [?1  -  ?A]  -  BF  (cA  -  C5) 

+  R(t)  [cA  -  C2] 

+  (CA  “  al>  +  f2  '^A  -  V  +  ^  W 

Then,  the  equation  for  linear  motion  of  the  recoiling  parts  is  (by 
substituting  for  in  Equation  86  ) 

Equation  91  (See  Page  64) 

Equations  89,  90,  and  91  are  the  governing  motion  equations  which  define 
the  variables 

ip  .  fi  and  n 

(and  their  derivatives)  during  the  firing  cycle. 
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EQUATION  83 


H[Yf+Z*+  i\A  +  +  2(nAYT+  iAZT)cos(j+B)  +  2(nAZT-tAYT)siT>(j  +  6)] 

'<$  '  +  MelY?4Z*  +  n*+(*+  +  LbZt)C0S(j  +  6)+  2UlsZr~  tBYT)s;ner  +  6)) 

+  +IIA  +  4^4] 

Ma  +  <1  +  <  V'r  +  Ijfrf-ostr  +6)  +  (r,AZT  -  iAYT)slnir  +  0)] 

+  $  •  +Me[ril  +£  4-  (ngYr-t  i^^oslj  +  6)  -+-  l\2T-  tjjslnlj  +  B)} 

+  Ua  +  Ie1  j 

•  *  f  \ 

+  ’M-MJC-  YTsin(j+6)  +  ZTeos<r  -fflJjj  - 
-  Yrsin( 7  +  e>  -h  Zrcos( 7  +  e>]  -  cA> 


"aU+  (H+  inAYT+iA2, 
+  1  Ia+V 
+  e  {»aWa  +  il)  + 


+  IaI-MaW  “ 


[6(fj]  i,  -  e  0  -  *A[r)Acos<T  +  6 


-  WA  [  rrcos  4>  —  ZTSin  <t>  t  iAcos(  7  +  6  -E  <t»  —  tA  s< »  (7  f  8  +  $)] 


-  WB[YTcos<t>  -  ZTsin$  +  if  cos<7+  6  +  <b)  —  Cpsir>  1  7  ->r  6  +<t>)] 


—  lv0[K0cos<t  -z0siv<t>]  —  K6L  <  <t>  -  Vstat'c  > 


+  2&B  I  ^[fn-yr  +  lAZr)sm(y +B)-(nAZT- t.AYr)cos<7  +  6)} 


[+  «8[fVr  +  (B^s‘n(Y  +  B)  -  <nBZT  ~  ieYrl  COSI7  1  B>] 


2BVa  Ma  [  vA  +  Yr  coslj  +  6)  +  ZTsin( 7  +  0)} 

—  2  qA i>  Ma [  nA  +  Ytcos(7  -r  A)  +  Z7  sin(7-r  8>) 


.2 1  w,[ru4rT  +  ii:zT)sin (j  +  e<  -  mAzT-  iayt> cosit  +  *>) 
|-f  Mg {( HB Yt  +  iBZr)sin(7  +  9)  -  UfBZT  -  {/TJCOSf7  +  9!) 


+  2  K; 


[EU  t  2£L2cos(y-i  91 
-[ELI  -*-2 


I  L—  [ctj - 


I  +  2EL2  cosIt  +  9) 
EU  -r  2EL2COS  7  - 


■f- 4>e[o]  +  «%!-2Wj 


2  j  Ma[  1^  + 


+HBi<%yT  +  ' 


R(t)-  BU)-[H'\S,\+  U^SJ+Fjsgr.l^)  -  BE  -  WAsin(j+  9+<t»  +  2Q6MAnA  +  6*MAn 


EQUATION  90 


\J4  8)  4  2 (rtAZT  -  iAyT  j  shur  +  e)] 
(J4-0)  +  Zfri BZT~  iBYT)sm<j-h8)] 


tf+6)  +  (\ZT-iAYT)stmy  4  Si] 
fVr~  tBYr)si»(7  +  8)) 


Ej-co sty  +  Wj|  - 


Ma[t^4-  CA  +  (riAYr  +  iAZr)cos<j  4  SI  4  (nAZr  -  iAYT)siniy  +  fij] 

£  ■*-  uBlnl  +  cb  +  (naYr  +  iBzTicos(r  +  ei-h cvBzT  -  teYT>sin  n  +  ei) 

■M4+V 


+  «*  {*(.[’! 1  +  O  +  +  0  +  t/4  ■#■/,]} 


[flfN]  (t  -  c.0  -  WAlnAcostr  4  8  4<t>)  -  ( Asin(j  4»4<t»] 


I  -+0)  —  C  sin  (j  -h  e  +$)] 


7+8  +-<t»-tasiTU7  -he  4(f»] 

o 


-  :AYT)cos(y  + 

-  ieYT)  cos( 7  4  6)]j 

sin  (7  4-  8)  j 

•h  ZT  sinij  4  6)} 

4-8)-  <nAZT-  iAYT)C0S(j  4  SJ] 
4-8)—  (HgZr  —  lBYT)  cosf7  4  8J] 


—  HJ, [nec°sf7  4  8+  <t»  -  CgSinir  4  8  4<t»] 


j" [ELI  4  2EL2cos(,  +8)  -  2EI.3 sinij 


-[EI.I+2EL2]  j  “j 

( _ EL2 sin  (7  +6)  +  EL3cos(y  +  6) 

[EU+2EL2cos(r-h  0)  -  2EL3sin(j  4  8)]1'* 


+  2{K 


[EL I  4  2EL2 cos( 7  +  8)  -  2EL3 sin (7  4  8)]"*' 
+  2EL2ccs  7  —  2EL3  sin  7] ui  -  r 


EL2sm  (74-8)  4  EL3cos(j+9) 
[EU+2EL2cos (7  +8)-  2EL3sm(7  4  $)]"* 


+  te[o]  +e%[-2MA  T(j  +  *J>l-2MAnA] 


•z  f  Ml[iriAYT+  iAZT)sin(y  +8)  —  (*AZT  -  (4>'rJcosf7  +  S)]] 
\+  Mb[(%Yt  +  {gZJSin<7  4-6)  —  tleZT  -  LYr)C0S(7+  gj]j 


EQUATION  91 


Bit)  —  [m(IS,I  +  WjpISgi  +  F6)sgn(jiAi  ~  8F  ~  WAsin(7  4  8+  4>)  4  2$8MAnA  4  82MAr\A  4  4?MA[riA4  YTC0S<7  4  8)  4  ZTsin(y  49)} 


O 

C^' 


Note  that  the  breech  force  [B(t)]  is  the  only  forcing  function  | 
included  in  Equations  39  and  90.  However,  the  following  functions  , 
are  included  in  Equation  91:  j 

i 


1.  R(t)  -  Force  due  to  recoil  mechanism 

2.  -  [f^  +  f2  +  Fg)  sgn  (n^)  -  Force  due  to  guide  friction 

3.  BF  -  Force  due  to  forward  buffer 

These  must  be  defined  directly  or  in  terms  of  the  motions 
which,  produce  the  forces. 

Previously,  the  force  due  to  the  recoil  mechanism  [R(t)]  has 
been  defined  (Equation  19)  in  terms  of  the  oil  pressures  ,  P2 
and  P^  (Figure  2) .  A  comparison  of  the  variables  n  and  x  leads 
to  the  following  definitions 


(92)  x  -  nAo  -  nA 


(93)  i--n. 

H 


,  .  •«  •• 

(9->)  x  »  -  n  . 


These  relations  can  be  used  to  solve  the  equations  (Table  I)  of 
fluid  flow  in  the  recoil  mechanism  simultaneously  with  the  equations 
(Equations  G9,  90,  and  91)  defining  carriage  response  to  firing  loads. 

In  a  similar  manner,  the  effect  of  the  forward  buffer  (required 
in  case  of  a  misfire)  can  be  included.  The  necessary  equations 
are  derived  in  the  Appendix. 

Ir.  either  Equation  17.2  or  Equation  21,  the  guide 
was  assumed  to  be  constant  and  denoted  by  the  term  F  . 
accurate,  this  force  is  dependent  on  the  reactions  (Si 
applied  at  the  guide  bearing  points  (Figure  12).  That 
the  carriage  response  model,  guide  friction  is 

(fl  +  f2  +  V  8S»  (x)“  ~  I^i !  sx  i  +  u2|S2!  +  Fg]  agn  nA 

The  reactions  and  are  determined  in  the  following  section. 


friction 
To  be  more 
and  S2) 
is,  for 


si-jjJ 
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EVALUATION  0?  DYNAMIC  REACTIONS 


In  reality,  this  is  a  continuation  of  the  carriage  response  model. 
In  fact,  the  clip  reactions  (S^  and  S;>)  wist  be  determined  to  solvci 
Equations  S9,  90  and  91.  These  are  evaluated  as  follows: 

From  Equation  87 


S5l  +  S2  =  WA  cos  (y  +  6  +  4>)  +  Ma  Q 

»  • 

or,  substituting  for 

S-  +  s.  -  V)  cos  (y  +  0  +  <J>) 

1  2  A 

+  K  Jv  {Yt  cos  (y  +  0)  +  ZT  sin  (y  +  8)  + 

+  0  -  2  q>6  5^  +  2  6  fi^  +2  $ 

-  02CA  -  i2[cA  -  Yt  sin  (y  +  9)  +  ZT  cos  (y  +  0) ] J 

I 

From  Equation  88  , 

1 

S1  ^  -  52  (d  -  a)  -  B(t)  [cx  -  <A]  -  SF  (tA  -  C$) 

+  R(t)  [cA  -  c23 

^A  “  ai^  +  f2  “  °2^  +  “  aG'  *  G^  ' 

-  -  I 

I 

-  IA  (6  +  ♦)  ' 

i 
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(95)  RAIL  -  -  [m1IsiI(Ca  -  +  u2|S  2I<CA  -  a2)  +  (cA  "  FG]  sgn  (  „ 


(96)  XYZ1  -  B(t)  [C1  -  CA]  ~  BP  (^  -  t5) 


+  R(t)  U.  -r  ]  -  I  (V  +  V) 

A  Z  A 


’4>  [Y^  cos  (y  +  9)  +  ZT  sin  (y  +  8)  +  nA] 

(97)  XYZ  -  M.<  +  *8n  -  2$k  +  28n  +  2n  A  -  8  2  cA 
A  A  A  A  A  A 


-  <}>  (C  -  Y  3 in  (y  +  8)  +  Z  cos  (y  +  0)] 
A  I  1 


+  WA  cos  (y  +  8  +  $) 


Then, 


S,  +  S„  -  XYZ 
1  2 


Sl(£)  -  S2  (F  -  a)  -  XYZ1  +  RAIL 


Solving  simultaneously, 


(98)  S, 


XYZ  (a)  -  XYZ1  -  RAIL 


(99)  S  -  XYZ  -  S2 


These  may  be  solved  for  Sj_  and  S2  by  a  simple  iterative  scheme  after 
assuming  starting  values  and  continuing  until  the  solution  converges. 
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•  c: 


The  force  in  the  two  equilibrators  has  been  defined  *•  (page  55) 


(76)  my  +  0-  2[PEQ  -  KEQ  (LEL,  +  Q  -  LEL^J 

and  the  force  in  the  two  elevation  struts  is  given  by  (See  Figure  15) 

EEL.  *  2£  (A  -  r)  K  +  P  ] 

Y  +  6  s 

or 

(100)  FEL  “  2[(LEL  ,  „  -  LELy  -  r)  K  +  P  ] 

y  +  0  Y  +  9  3 

where  (page  53) 

(75)  LEL  ,  -  [ELI  +  2EL2  cos  (Y  +  9)  -  2EL3  sin  (Y  +  &))±/2 

Y  +  o 

In  evaluating  the  trunnion  reactions,  the  components  of  acceleration 

for  the  mass  M  are: 

B 

in  the  n  direction, 

R  «  B  cos  (y  8  +  $)  +  C  sin  (y  +  8  +  $) 
nB  P. 

and  in  the  C  direction, 

R,  -  -  *B  sin  (Y  +  9  +  $)  +  c  C08  (Y  +  9  +  4>) 

**  B  B 


68 


Substituting  for  BB  and  C3  and  simplifying  allows  writing  these  components 
in  the  following  form 

Rn  -  -  V  [SB  -  Yx  sin  <Y  +  0)  +  ZT  cos  (Y  +  9)]  -  8*CB  -  2$% 

•  2  *2 

i  -  $  l^g  +  Yt  coa  (y  +  0)  +  Zj  sin  (y  +  0)]  -  9  ig 

and 

R  -  V  IhB  +  Yt  cos  (Y  +  0)  +  ZT  sin  (Y  +  0)J  +  0*nB  -  2$0CB 

•  2  *2 

-  4  [tB  -  Yt  sin  (Y  +  9)  +  ^  cos  (Y  +  B)]  -  9  Cg 

Then,  from  the  free  body  diagram  of  the  cradle  [Figure  12(b)] 

lFn  -  -  Th  -  FEQy  +  e  cos  t  -  R(t)  +  BF 

+  FEI^  +  g  cos  T  -  WB  sin  (Y  +  9  +  ♦) 

+  +  f2  +  Pol¬ 

and 

^  -  HgRj.  -  Tv  +  FEQy  +  0  8in  T  “  S1  "  S2 

-  FELy  +  0  sin  t  -  Wg  cos  (Y  +  9  +  * ) 


where  (Figure  12) 


Z  —  7 

-1  1  EL  .  ,  . 

i  ■  ten  — - k—  +  (y  +  8) 

tEL  "  1 


NOTE:  The  angle  t  will  be  in  the  first  or  second  quadrant 


,  Z  +  n_.  sin  (y  +  8)  +  ?  cos  <y  +  9)  -  Z„ 

<101)  '  *  “•  +  •> 


Then, 


(102)  TH  -  -  Mg  I  V  [CB  -  Yt  sin  (y  +  6)  +  ZT  cos  (Y  +  8)]  +  8cB  +  2*6^ 


+  $  InB  +  008  (y  +  8)  +■  ZT  sin  (y  +  8)]  + 


•s] 


+  FEQY  +  0  C08  T  +  R<c>  -  U1lsif+  »2,S21+  FCJ  88n  ^ 


-  BF  -  PEL 


(103)  T 


Y  +  6  cos  T  +  Wg  sin  (Y  +  8  +  ») 


Mg  ^*tnB  +  ^  cos  (y  +  6)  +  ZT  sin  (y  +  0)]  +  VnB  -  2£ecB 


-  $  ~  sin  <Y  +  8)  +  ZT  cos  (y  +  8)]  - 


~  *®QY  +  0  8in  T  +  Si  +  s2  +  FEL^  +  0  sin  t  +  Wg  cos  (y  +  8  +  <fr) 
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The  free  body  diagram  of  the  mass  is  shown  in  Figure  16.  Since  this 
mass  is  also  assumed  to  simply  rotate  around  point  0,  the  constants. 


(104) 


S  sltan 


-i  S 


(105)  Ru  .  [Yc2  +  2d2]' 


may  be  used  in  evaluating  the  reactions  at  the  support  points.  As 
defined  on  Page  5<1  the  total  spring  deflection  is  given  by 


AG  "  %  +  AG 


static 


Then  ,  since 


%  “  '  L* 


‘  static 


static  r . 


(106)  » 


static 


static 
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Where  P8tatic  be  evaluated  as  (Page  50) 


ba  +  WA  +  bb  wb  *  bd  wd 


static 


Therefore, 


ag  "  -  L  (*  -  *statlc> 


The  reaction  at  the  effective  spring  is  given  by 


Vspring  "  KG  AG 


( 10?)  VsFting*"KG  L  <*  "  Wic> 


the  free  body  diagraa  or  the  Hass  Mp  {Figure  16) 


EFg  -  sin  (6  +  <P) 


Tfi  cos  (y  +  9  +  ♦)  -  Tv  sin  (Y  +  8  +  *) 


FTXy  +  0  cos  (t  -  y  -  6  -  <t>) 


’  *  9 


(t-y~  8  -  $)  “  H 


.. 


EFc "  %  M  cos  ^ +  ^  " 


V  J  +  FEL  „  sin  (t  -  Y  “  9  -  <>) 

spring  y  +  9 


-  rH  sin  (\  +  e  +  $)  -  tv  cos  (y  +  a  +  $) 


-  FEQ  sin  (t  -  Y  “  9  -  $) 

Y  +  9 


-  W  +  V 

D  pivot 


Therefore, 

(108)  V  .  «  Tu  sin  (y  +  9  +  *)  +  T„  coa  (y  +  0  +  *) 

pivot  H  V 

-FEL  4  ain  (t  -  y  -  0  -  ♦) 

y  4-  6 

+  ^  +  9  sin  (t  -  y  >  9  -  ♦>  +  WD  -  V8prlng 

+  &d  C0B  (B  +  V 

(109)  H  -  TH  COB  (Y  +  6  +  ♦)  -  Ty  sin  (y  +  8  +  $) 

4-  FEL  cos  (t  -  y  “  6  ~  *) 

Y  4-  6 

-  FEQ^  +  Q  cos  (t  -  y  -  9  -  ♦)  -  [M^  Rjj  sin  (F  +  $)] 
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EVALUATION  OF  STATIC  REACTIONS 

Under  static  conditions,  the  forces  B(t) ,  BF,  and  f2  shown 
in  Figure  12a  will  be  assumed  equal  to  zero.  Then, 

£Fn  -  0  -  R(t)  -  WA  sin  (Y) 

R(t)  -  WA  sin  (Y) 

-  0  -  S1  +  s2  -  WA  cos  (Y) 

+  $2  *  WA  cos  (Y) 

“:<Kass  Center"  0  =  Si  (a>  -  S2  tt>  -  1}  -  R(t)  I?A  -  C2J 
Sx  (a)  ■*  S2  (b  ■“  a)  *  WA  [CA  ~  T 

Solving  simultaneously 

W  [(a)  cos  Y  -  (C.  -  C7)  Sin  y] 

(no)  s2  -  — - r — - 

b 

(111)  “  WA  cos  Y  "  s2 

The  force  in  both  equilibrators  is  determined  in  the  manner  described  in 
developing  the  basic  model.  Therefore, 

(112)  FEQ  »  2  {PEQ  -  KEQ  (LEI^  -  LEL^)  ] 
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To  define  the  static  loud  In  the  elevating  mechanism  (wh^ch  was  considered 
approximately  zero  when  developing  Models  A  and  B),  on  examination  of 
Figures  11  and  12  shows  that,  for  the  tipping  parts  (M^  and  Mg) 

£M  *  0  »  W  (n  cos  Y  —  C  sin  Y) 

T  a  AO  A 

+  WD  (n  cos  Y  -  4  sin  Y) 

B  B  B 

-  FEQy  cos  T  [Cel]  -  FEQy  sin  t  [n^] 

+  FELy  sin  T  [na]  +  FELy  cos  t 

Under  static  conditions, 9  and  $  are  equal  to  zero.  Therefore, 


(113)  t  -  tan 


ZT  +  "EL  Sln  Y  +  CEL  COa  Y  ~  ZEL 

YEL  -  V  nEL  C0S  Y  +  SL  8ln  Y 


and  the  loaU  in  the  elevating  strut  is  given  by 


Ui4) 


FEL 

Y 


FEQy  (nEL  Sin  T  +  CEL  COS  T)  '  WA  (  °AC  C0S  Y  "  CA  9in  Y> 

**  ^g  (rig  cos  Y  "Cg  Y^ 

n£L  sin  t  +  ;EL  cos  “ 
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In  a  similar  manner,  Che  following  equations  allow  for  determination 
of  the  trunnion  ’■aactiona. 


IF 


7U  -  FEQ  cos  t  +  FELy  cos  t 


-  W  sin  y  -  W  sin  y 
A  B 


IF. 


T  +  FEQ  sin  x  -  FEL  sin  x 

V  Y 


-  cos  y  “  W  cos  y 

A  B 


Therefore, 


(115)  T  -  (FEQ  -  FEL  )  cos  x  +  (W^  +  W  )  sin  y 

H  v  Y  *  “ 

(116)  Tv  -  (W  +  W3)  cos  y  +  (FELy  -  FEQ^)  sin  t 
Finally,  from  Figure  11 


EM  -  0  -  V  L  -  W.  [Yt  +n cos  Y  -  ?A  sin  y] 

U  spring  A  T  AO  A 


WB  [Yt  +  nB  cos  Y  -  CB  sin  y] 


UD  [V 


IFh  «  o  •  h 
IF..  - 


! 


I 


■  I 

1  j 


0  -  V 


W.  -  W  -  W.,  +  V 


Therefore,  under  static  conditions 


(117)  H  -  0 


Wa[Yt  +  nAQ  cos  y  -  £A  sin  yl  +  Wg[Y^  +  nB  cos  Y  “  CB 


+  W  Y 
D  D 


(118)  V 


spring 


(119)  V  .  „  =  W.  +  W„  +  W_  -  V  , 

pivot  A  B  D  spring 


When  the  recoiling  mass  returns  to  che  latch  position  at  the  end 
or  any  firing  cycle  or  when  It  reaches  the  limit  of  forward 
travel  in  the  case  of  a  misfire,  additional  relative  motion  between 
the  masses  MA  and  Mg  is  prevented  by  mechanical  stops.  Consequently, 
the  preceding  system  of  equations  must  be  modified  by  eliminating 
Equation  91  completely  and  holding  the  value  of  0  (and  x) 
constant  from  this  time  on.  At  the  same  time,  in  the  remaining 
equations , 


sgn  (nA)  -  0 


sin  y] 


78 


REFERENCES 


1.  CARVER,  Harvey,  "Final  Report  -  An  Investigation  of  the  Firing  Out 
Of  Battery  Principle  on  a  Field  Artillery  Weapon",  Report  64-1348, 
Research  &  Engineering  Directorate,  Rock.  Island  Arsenal,  April  1964, 

2.  TOWNSEND,  Philip,  "Analog  Computer  Study  of  the  Initial  Momentum 
Technique  in  Reducing  Rod  Pull  Requirements  for  a  105na  Howitzer" 

Report  64-2615,  Research  &  Engineering  Directorate,  Rock  Island 
Arsenal,  September  1964. 

3.  PETERSEN,  Donald  F. ,  "Exploratory  Development  Studies  of  105nm 
Fire  Out  Of  Battery  Weapon",  Report  68-799,  Research  &  Engineering 
Division,  Pock  Island  Arsenal,  March  1968. 

4.  SEAMANDS ,  Robert  E. ,  "Exploratory  Development  of  Howitzer,  Light, 

Towed;  105mm  Soft  Recoil,  XM204" ,  Artillery  Systems  Lab  Report 

RE  TR  70-179,  Artillery  Systems  Laboratory,  Research  &  Engineering 
Directorate,  U.  S.  Army  Weapons  Coiamand,  Rock  Island  Arsenal, 

July  1970. 

5.  WHITCRAFT,  James  S.,  "Engineering  Design  Test;  of  Howitzer,  Light, 

Towed:  105mm  Soft  Recoil  (Safety  Certification)",  Aberdeen 
Proving  Ground  Report  APG-MT-3664,  October  1970. 

6.  JOHNSTON,  Major  Franklin,  Jr.,  "Military  Potential  Test  of  Howitzer, 
Light,  Towed:  105mm  Soft  Recoil,  XM204",  U.  S.  Army  Field  Artillery 
Board,  Final  Report,  28  May  1971. 

7.  NERDAHL,  Michael  C.  and  FHAN1'?,  Jerry  W. ,  "Prediction  of  System 
Motion  Based  on  a  Simplified  Mathematical  Model  for  a  Soft  Recoil 
(firing  Out  Of  Battery)  Mechanism",  Technical  Note  Art  3-69, 

Research  &  Engineering  Directorate,  U.  S.  Army  Weapons  Command,  May  1969. 

8.  NERDAHL,  Michael  C.  and  FRANTZ,  Jerry  W. ,  "Engineering  Analysis  XM46 
Recoil  Mechanism,  Design  of  Functional  Controls  and  Prediction  of  System 
Motion",  Artillery  Systems  Lab  Report  69-1398  w/oupplement,  June  1969, 
Artillery  Systems  Laboratory,  Research  &  Engineering  Directorate, 

U.  S.  Army  Weapons  Command. 

9.  NERDAHL,  Michael  C.  and  FRANTZ,  Jerry  W. ,  "Engineering  Analysis, 

XM204  Howitzer  (Soft  Recoil) ,  Effect  of  Carriage  Flexibility  on 
Weapon  Response",  Artillery  Systems  Lab  Report  RE  69-2779,  Artillery 
Systems  Laboratory,  Research  &  Engineering  Directorate,  U.  S.  Army 
Weapons  Command,  December  1969. 

10.  GRITT0N,  Robert,  "Initial  Engineering  Design  and  Test  Effort  on 
the  XH204  Advanced  Development  Model",  Technical  Note  A  2-71, 

Artillery  and  Air  Defense  Weapons  Systems  Directorate,  Rock  Island 
Weapons  Laboratory,  U.  S.  Army  Weapons  Command,  August  1971. 


79 


SYMBOL  TABLE 


While  symbols  and  subscripts  have  been  defined  (formally  and/or 
pictorially)  as  they  were  introduced,  the  following  list  is  included 
as  a  convenient  reference.  Model  complexity,  number  of  models, 
and  use  of  standard  notation  has  resulted  in  some  symbol  duplication. 

'his  table  may  aid  in  resolving  questions  resulting  from  such  duplcation. 
in  all  models,  Che  basic  unite  are: 

inches 

pounds 

seconds 

radians 


A 

Subscript  showing  relation  to  mass 

*1 

Area  at  pressure  P^ 

Fig.  2 

A2 

Area  of  spear  buffer  -  pressure 

Figure  2 

A3 

Area  at  pressure  P^ 

Figure  2 

A4 

Effective  area  of  floating  piuton  (oil  side) 

Fig.  2 

N 

Effective  area  of  floating  piston  (gas  side) 

Fig.  2 

Area  of  each  recoil  rod 

Fig..  2 

*b 

Effective  area  of  each  front  buffer  piston 

Fig.  A-l 

A(t) 

Summation  of  forces  causing  positive 
accelerations 

Fig.  7 

al 

Orifice  area  between  pressures  P^  and  P^ 

Fig.  2 

*2 

Orifice  area  between  pressures  and  P ^ 

Fig.  2 

a3 

Orifice  area  through  velocity  sensor. 

Effective  area  between  pressures  and 

Fig.  2 

a3bf 

Value  of  a  before  firing 

*3af 

Value  of  a ^  after  firing 

®v 

Flow  area  through  spear  buffer  check,  valve 

Fig.  2  6 
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a  Flow  area  due  to  clearance  between  Fig.  5 

ea  spear  buffer  and  piston  head 

a(x)  Variable  orifice  area  dependent  on  position  Fig.  3 

of  spear  buffer 

a  Orifice  area  for  front  buffer  Fig.  A-l 

B 

X  Locates  centroid  of  area  under  B(t)  Figs.  8  &  9 

"a  Distance  from  rear  clip  reaction  to  mass 

center  of  recoiling  parts  Fig.  12  (a) 

B  Subscript  showing  relation  to  mass  Mg 

B  Measurement  parallel  to  0  -  B  axis  Fig.  11 

BF  Retarding  force  of  front  buffer 

B(t)  Breech  force 

b  Distance  between  front  and  rear  Fig.  12  (a) 

clip  reactions 

C  Measurement  parallel  to  0  -  C  axis  Fig.  11 

c^  Damping  coefficient  related  to  $ 

cQ  Damping  coefficinet  related  to  9 

Discharge  coefficient  for  orifice 

D  Subscript  showing  relation  to  mass  M  ^ 

D(t)  Summation  of  forces  causing  negative  Fig.  1 

acceleration 

Value  of  D(t)  at  t  »  0  Fig.  8 

D^  Value  of  D(t)  at  t  *  t  Fig.  8 

Dr  Constant  level  for  D(t)  during  recoil  Fig.  8 

Dcr  Maximum  value  of  D(t)  during  counterrecoil  Fig.  10 
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EQ 

Subscript  relating  variable  to  equilibrator 

EL 

Subscript  relating  variable  tc  elevating  struts 

F* 

Generalized  force  (torque)  causing  rotation  $ 

Fe 

Generalized  force  (torque)  causing  rotation  6 

FEL 

Force  on  elevating  struts 

Fig.  12 

FEQ 

Force  on  equilibrators 

Fig.  12 

NOTE:  Elevating  strut  and  equilibrator 

are  considered  to  be  separate  but  parallel 
units 

Ffp 

Packing  friction  for  floating  pi3ton 

Fig.  2 

F 

P 

Packing  friction  for  recoil  piston 

Fig.  2 

fg 

That  portion  of  guide  friction  which  is 
independent  of  clip  reactions 

Fig.  2 

fl 

Guide  friction  due  to  clip  reaction  S^ 

Fig.  12 

f2 

Guide  friction  due  to  clip  reaction  S2 

Fig.  12 

G 

Subscript  relating  variable  to  ground  spring 
or  recoil  guides 

9 

Acceleration  due  to  gravity 

S(v£) 

Pressure  drop  across  '  i'th  orifice 

Eq.  1 

H 

Step  function 

Eq.  20 

H 

Horizontal  reaction  at  ground  pivot 

Fig.  11 

U 

Constant  related  to  elevating  strut  load 

Fig.  15 

I  j 

2 


I 


■a 

S 

I.D. 

K 

K 

EQ 

K3 

“g 

k 


Impulse  -  Area  under  B(t) 

Mass  moment  of  Inertia  for 

Mass  moment  of  inertia  for  Mg 

Mass  moment  of  inertia  for  M^ 

Duration  of  ignition  delay 

Spring  rate  for  elevating  struts  Fig.  15 

Spring  rate  for  equilibrator  Fig.  13 

Spring  rate  of  each  front  buffer  return  spring  Fig,  A-l 

Effective  spring  rate  of  ground  spring  Fig.  11 

Gas  constant  -  Ratio  of  specific  heats 
for  gas  at  pressure  PN 


L  ' Y '  coordinate  of  ground  spring 

LELq  Length  of  elevating  struts  at  elevation  -  0° 

LEL  Length  of  elevating  struts  at  initial 

Y  elevation  (y ) 

LEL  „  Length  of  elevating  struts  at  rotation  6 

y  +  0 


*P 

H 

e; 

ma 


Mass  of  recoiling  parts  without  floating  Fig.  2 

piston 

Mass  of  floating  piston  Fig.  2 

Effective  mass  of  recoiling  parts  Eq.  27 

Mass  of  recoiling  parts  used  in  three-degree-  Fig.  11,  12 
of-freedom  model 

Mass  of  elevating  (but  non-recoiling)  portion  Figs  11  &  12 
of  the  weapon 

Maas  of  non-elevating  Portion  of  the  weapon  Fig*  H  &  16 
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N  Nissber  of  recoil  cylinders 


0 

Ground  pivot  -  Origin  of  0  -  BC  and  0  - 
coordinate  systems 

YZ 

Fig.  11 

P1 

Oil  pressure  in  recoil  cylinder  -  spear 
chamber 

buffer 

Fig.  2 

4 

P2 

Oil  pressure  in  recoil  rod 

Fig.  2 

P3 

Oil  pressure  in  recoil  cylinder 

Fig.  2 

P4 

Oil  pressure  in  recuperator 

Fig.  2 

PN 

Gas  pressure  in  recuperator 

Fig.  2 

P 

o 

Initial  gas  pressure  in  recuperator 

Fig.  2 

PR 

Oil  pressure  in  front  buffer 

Fig.  A-l 

PEQ  Preload  in  equilibrator  spring  (i.e.  -  at  y  «  Q*) 

~Q  Vector  locating  mass  center  of  recoiling  parts 

\  Vector  locating  mass  center  of  elevating  (but 

non-recoiling)  parts 

R  Rod  pull  -  Force  on  recoil  rod 
R(t) 

RD  Distance  from  ground  pivot  0  to  mass  center  of 
non-elevating  parts 

r  Constant  relating  tn  elevating  strut  load  Fig.  15 

S,  Normal  reaction  at  rear  support  of  recoiling  parts  Fig.  12 

Normal  reaction  at  front  support  of  recoiling  parts  Fig.  12 
sgn(u)  Algebraic  sign  of  the  variable  u 
S/  Preload  in  each  front  buffer  return  spring 
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T 

Elevation  trunnion  -  origin  of  T  -  n  C 
coordinate  system 

Fig.  U 

T 

Kinetic  energy 

T 

Duration  of  counterrecoil  control  function 

Fig.  10 

th 

Trunnion  reaction  -  1  1  to  T  -  n  axis 

Fig.  12 

Tv 

Trunnion  reaction  -  1 1  to  T  -  ;  axis 

Fig.  12 

t 

Time  variable 

t 

r 

Time  of  recoil 

Cf 

Time  of  firing 

V 

Potential  energy 

y 

pivot 

Vertical  reaction  at  pivot  point  0 

Fig.  11 

Vspring 

Vertical  reaction  at  ground  spring 

Fig.  11 

Recuperator  gas  volume  for  displacement  'x' 
Initial  recuperator  gas  volume  (x  ■  0) 

Fluid  velocity  through  '  i'th  orifice 


U 


eff 


Weight  of  reciling  parts  (without  floating 
piston) 

Weight  of  floating  piston 

Effective  weight  of  recoiling  parts 

Weight  of  recoiling  parts  for  three-degree- 
of-freedom  model 

Weight  of  elevating  (but  non-recoiling) 
portion  of  weapon 

Weight  of  non-elevating  parts 


Eq.  28 
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x,x,x  Displacement,  velocity  and  acceleration  of  Pig.  2 

recoiling  parts 

x  Recoil  displacement  for  which  spear 

s"  buffei  becomes  effective 

x  Recoil  displacement  at  which  spear  buffer 

v  check  valve  ceases  to  be  effective 

x  Recoil  displacement  at  which  counterrecoil 

n  control  by  the  spear  buffer  begins 

x  Recoil  displacement  at  which  front  buffer 

is  actuated 

x  Constant  counterrecoil  velocity  governed 

mir‘  by  flow  through 

x^  Firing  velocity 

x  Recoil  displacement  at  which  terminal 

'  counterrecoil  velocity  is  reached 

Terminal  counterrecoil  velocity 


Y  Measurement  parallel  to  0  -  Y  axis  Fig,  11 

y,y,y  Displacement,  velocity  and  acceleration  of  Fig.  2 

floating  piston 

Z  Measurement  parallel  to  0  -  Z  axis  Fig,  11 

a  t  coordinate  of  Fig.  12 

6  C  coordinate  of  f^  Fig.  12 

Angle  locating  mass  center  of  non-elevating  Fig.  16 

parts 

Y  Initial  angle  of  elevation  Fig.  11 


A  Deflection  of  ring  springs  in  elevating  struts  Fig.  15 
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C i 

r 

Specified  rise  and  fall  Claes  for  D(t) 

Fig.  S 

C 

Measurement  parallel  to  T  -  C  axis 

Pig.  11 

n 

Measurement  parallel  to  T  -  i  axis 

Pig.  11 

VVnA 

Defines  position,  velocity,  and  acceleration 
of  mass  M  with  respect  to  the.  elevation 
trunnion  A 

Pig.  12 

nA0 

Initial  value  of 

•  •• 

0,6, e 

Defines  relative  rotation,  velocity 

and  acceleration  between  M„  and  M  around  the 

a  D 

elevation  trunnion 

Fig.  11 

u 

Coefficient  of  friction 

O’ 

Density  of  recoil  oil 

T 

Angle  between  elevation  strut  and  the 

T  -  n  axis 

Pig.  12 

•  •• 

$  » 4> 

Defines  rotation,  velocity  and  acceleration 
of  mass  M  around  an  axis  of  rotation  fixed 

In  the  ground  at  point  0 

Fig.  11 
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APPENDIX 

MATHEMATICAL  MODEL  OF  FRONT  BUFFER 
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In  case  of  a  misfire  or  accidental  tripping  of  the  latch 
mechanism,  the  recoiling  parts  will  be  driven  to  their  forward 
limit  of  travel.  Consequently,  a  hydraulic  buffer  is 
incorporated  in  the  design  to  bring  the  moving  mass  to  rest  and 
protect  the  weapon  against  impact  loading.  Since  fluid  flow 
from  the  recuperator  is  restricted  after  the  velocity  sensor  has 
functioned,  the  front  buffer  design  will  be  based  on  the 
maximum  expected  firing  velocity.  This  mechanism  is  shown 
schematically  in  Figure  A-l. 

To  limit  the  initial  impact  loading  which  occurs  as  the  recoiling 
mass  (Mj^)  hits  the  buffer  pistons  (Mg)  ,  a  spring  assembly 
(rate  *  K^)  was  added  to  the  recoiling  parts.  The  spring  rate 
was  chosen  so  that  the  velocities  of  Mr  and  Mg  are  approximately 
equal  when  the  spring  deflection  becomes  equal  to  £ 

After  this  time,  the  resisting  force  of  the  front  buffer  is 
defined  by: 


BF  -  2AbPb  -  2  K3(x  -  Xg)  +  2S 


where 


S*  ■  Preload  in  each  buffer  spring 
Kg  “  Spring  rate  of  buffer  spring 

=*  Recoil  displacement  for  actuation  of  front  buffer 

u 

Ag  *  Effective  area  of  each  buffer  piston 


Pg  “  Buffer  oil  presasure. 


Now,  the  pressure  drop  across  the  orifice  aB  is  given  by 
2 

Cr  vb  ,  . 

®  *  Tg  ^7  3gn  (V 
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and,  since 


*b  *  “  "  Vb 

j  2  2 

.  ar  s,  . 

AP  “  -  Yg  — 2  x  Sgn  (x) 

*B 


For  a  more  complete  description  of  a  design  procedure  for  this 
mechanism,  see  Reference  8,  paces  97  -  117. 
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